Fluctuating Commutative Geometry 



Luiz C. de Albuquerque 1 , Jorge L. deLyra 2 , and Paulo Teotonio-Sobrinho 2 
(1 ) Faculdade de Tecnologia de Sao Paulo - DEG - CEETEPS - UNESP, 
Praqa Fernando Prestes, 30, 01124-060 Sao Paulo, SP, Brazil and 
(2) Universidade de Sao Paulo, Instituto de Fisica - DFMA 
Caixa Postal 66318, 05315-970, S&o Paulo, SP, Brazil 

We use the framework of noncommutative geometry to define a discrete model for fluctuating ge- 
ometry. Instead of considering ordinary geometry and its metric fluctuations, we consider generalized 
geometries where topology and dimension can also fluctuate. The model describes the geometry of 
spaces with a countable number n of points. The spectral principle of Connes and Chamseddine is 
used to define dynamics. We show that this simple model has two phases. The expectation value 
(n), the average number of points in the universe, is finite in one phase and diverges in the other. 
Moreover, the dimension 5 is a dynamical observable in our model, and plays the role of an order 
parameter. The computation of (S) is discussed and an upper bound is found, (8) < 2. We also 
address another discrete model defined on a fixed d = I dimension, where topology fluctuates. We 
comment on a possible spontaneous localization of topology. 

PACS numbers: 04.60.-m, 04.20.Cv, 04.50. +h, 02.10.Yn, 05.70.Fh 

Keywords: Quantum Gravity; Noncommutative Geometry; Random Matrix Theory; Phase Transitions. 



2 



I. INTRODUCTION 



A possible approach to quantum gravity is to consider it as an Euclidean quantum field theory. 
In the functional integral approach, one tries to compute a sum over metrics g^ defined on a d- 
dimcnsional manifold M. The dynamics is fixed by a choice of an action that is usually taken to 
be the Einstein- Hilbert action. An exact computation of such functional integral is not available 
and several approximated methods have been developed. The main idea is to use some kind of 
discretization in order to approximate the functional integral by a large number of ordinary integrals 
or sums. One such method is known as the dynamical triangulation and consist in summing over 
the triangulations of a given manifold. The situation becomes more complicated when one wishes 
to include quantum topology change. In this case, the manifold itself is also not fixed and a sum 
over topologically non equivalent manifolds M has to be included. The overall picture seems to be 
under control only for d = 2 where a classification of topologies is possible. Such a classification is 
proven to be impossible for d — 4 and it is an open problem for d = 3. 

The framework of noncommutative geometry is appropriate to explore some difficult questions 
in quantum gravity. We illustrate how noncommutative geometry can be used to generalize Eu- 
clidean quantum gravity, i.e, fluctuating geometry. Instead of considering ordinary geometry and its 
metric fluctuations, we consider the generalized geometries where besides the metric, topology and 
dimension can also fluctuate. 

The basic idea coming from noncommutative geometry is that one can describe a Riemannian 
manifold [M, Qpv) in a purely algebraic way. There is no loss of information if, instead of the data 
(M, g^), one is given a triple (A, H, D), where A is the C*-algebra C°(M) of smooth functions on 
M, Tt is the Hilbert space of L 2 -spinors on M, and T> is the Dirac operator acting on TC. From the 
Gelfand-Naimark theorem it is known that the topological space M can be reconstructed from the 
set A of irreducible representations of C°(M). Metric is also encoded, and the geodesic distance can 
be computed from T>. Here we will consider only commutative spectral triples - this is enough to go 
much beyond ordinary geometry. In particular one can treat all Hausdorff topological spaces in this 
way. Given a pair (M, <7 M „), one can promptly construct the corresponding triple (C°(M), L 2 (M), V). 
However, not all commutative spectral triples, or generalized geometries, come from a pair {M, g^ v ). 
Nevertheless one can always associate a Hausdorff space M — A to a commutative spectral triple, 
where A denotes the set of irreducible representations of A. However, the space M may not be 
a manifold. Once we trade the original Riemannian geometry for its corresponding commutative 
triple we need a replacement for the Einstein-Hilbert action Seh- The so-called spectral action of 
Chamseddine and Connes Q is one possible candidate. It depends only on the eigenvalues of T> 
(the spectral principle) and contains Seh as a dominant term. In this paper however we shall use 
another spectral action. 

The spectral action can be written for any triple, regardless of whether it comes from a manifold 
(M, g^ v ) or not. In the spectral geometry approach it is conceivable to write the partition function 



were the "sum" is over the set X of all possible commutative spectral triples and S depends on the 
spectrum of V. It includes all Hausdorff spaces and therefore all manifolds of all dimensions. 

Apparently, there is no advantage in considering the partition function since it is by no means 
easier to compute. However, the algebraic approach provides us with a natural way of defining 
discrete approximations for the theory. For that it is enough to replace the algebra A by a finite 
dimensional algebra A n . In this approach to discretization there is no need to introduce a lattice or 
simplicial decomposition of the underlying space. The approximation of A by a finite algebra works 




(1) 
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even if the spectral triple does not come from a manifold. In this sense, it gives us a generalization of 
ordinary discretizations |y]. The simplest discrete model one can consider is a simplification of ijTJl. 
where instead of summing over the set X of commutative spectral triples, we take a subset X C X . 
The set X consists of points x = (A, TL, D) where A is a commutative algebra with a countable 
spectrum A, i.e., with a countable number of irreducible representations. Therefore the underlining 
space M = A has a countable and possibly finite number of points. 

Most of the results discussed here were reported in In section 2 we reintroduce the discrete 
model of Q. In section 3 we show that this simple model has two phases. The expectation value 
of the number of points diverges in one phase and it is finite in the other phase. The definition of 
dimension in noncommutative geometry is recalled in section 4. A estimate of the expectation value 
of the dimension is discussed in section 5. In section 6 we briefly discuss another model for random 
geometry where the dimension is fixed and equal to one, but the topology fluctuates. Further results 
and details on this second model will be reported elsewhere |4|. 



II. DISCRETE MODEL 



The exact computation of (1) is a major goal, not yet accomplished. In this paper we discretize 
by sampling the set X with finite commutative spectral triples. We will think of it as a useful toy 
model, which seems to capture some of the main features of the full one, Eq. (1). For instance, the 
key role played by the eigenvalues of the Dirac (or Laplace) operator in the spectral action approach 
was emphasized in Q. In our model they are also the natural dynamical variables due to the 
connection with random matrix theory. An important ingredient of the model is that the number 
of points can fluctuate. Moreover, in our simple model the space-time dimension is a dynamical 
observable and its expectation value can be computable from first principles. 

Let us describe the ensemble X C X of geometries we will consider. A point of x £ X is a 
commutative spectral triple x = (A,H,D) where the commutative C*-algebra A has a countable 
spectrum A. We divide X into subspaces X n consisting of triples (A n ,H n , D) such that A n has a 
fixed number n of points. From the Gelfand-Naimark theorem it follows that elements of A n are the 
(possibly infinite) sequences a = (ai,a2, ...,a„), cij G C. The Hilbert space Ti n is given by vectors 
v = (vi, ...,v n ) with norm ||v|| 2 = Y17=i v i < 00 • The elements of A are represented by diagonal 
matrices a — diag(a!, ...,a n ) acting on Ti n . Finally, the operator D is a n x n self-adjoint matrix. 
We will sample the space X by Xi,X2, Xn and eventually take the limit N — > oo. 



III. DYNAMICS 



Let L be a length scale such that the operator T> given by T> = D/L will be the analogue of the 
Dirac operator. The Chamseddine-Connes action depends on a cutoff function of the eigenvalues of 
D/L. The cutoff function is zero for eigenvalues of D greater than L and one otherwise |2j,l?J. In 
other words, the Boltzmann weight in Eq.Q would be one outside a compact region in the eigenvalue 
space, leading to a divergent partition function (see Eq. H110 ). Let us consider a quadratic action 
instead: 

S[x]=Tr(?) 2 = (2) 
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where A is the inverse of Planck's length l p , and (3 — (l p /L) . Finally, we define the partition 



function Z N ({3) = J2n=o z n{@) where 

z n (/3) = / [dD]e-f™^ (3) 



is the partition function restricted to X n , in other words, an integral over all independent matrix 
elements Dij , where [dD] is the usual measure for n x n Hermitian matrices . The partition function 
z n ({3) defines the one-matrix Gaussian Unitary Ensemble 0]. A straightforward computation gives 

z n (P) = 2H^. (4) 

The expectation values of an observable O(Dij) restricted to X n and for the entire ensemble are 

(0) ntfl = J [dD]Oe-^ D ^/z n {(3) and (5) 

N 

(0)(/3) ee *£p(n,f3)(0) n ,p, (6) 



respectively, where the function 



P(n,(3) = ^%- (7) 



is interpreted as the probability of having a universe with n points. The simplest observable in our 
model is n, the number of points in A. By definition, n is constant in X n , therefore («) n ,/3 = n - 
Thus we get 

Wl/ 3 ) = — ~ ~^2-- ( 8 ) 

The mean (n) ("average number of points in the universe") is not a continuous function of (3 at 
j3 c = 7r/2, signaling the onset of a phase transition. Besides straightforward numerical calculation, 
there are other ways to show that the sum ||SJ| converges for (3 > (3 C and diverges for (3 < (3 C . 



IV. DEFINING DIMENSION S 

For (3 < (3 C the relevant universes have (n) — oo and An/{n) = 0. For a oo-dimensional D one 
can define the dimension 6 of the space A from the eigenvalues of D. Let {no(D), fii(D), ...} be the 
modules of the eigenvalues (i.e. the singular values) of D organized in an increasing order. By the 
Weyl formula jlj , the dimension S is related to the asymptotic behavior of the eigenvalues for large 
k: fik(D) w fc*. By definition 5 = for finite dimensional spectral triples. We can argue that (5) is 
of the form 



This follows from the fact that for (3 > j3 c the probability P(n, (3) is localized around some finite n. 
Hence (5) works as an order parameter. The value f3 c — ir/2 separates (S) — from the rest. 
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In order to study the dimension we need to consider the spectral (^-function 

n 

C(z)= lira J2^ Z = ^(\ D \~ Z l (10) 

k=0 

where D is an oo-dimensional matrix (/io > 0). The relation between the dimension and £(z) has 
been discussed in [ljj. For large enough values of a — Re(z), Tr(\D\~ z ) is well defined. One says 
that D has dimension spectrum Sd if a discrete subset Sd = {si, S2, ■■•} C C exists, such that £(z) 
can be holomorphically extended to C/Sd. This definition is consistent with the Weyl formula. The 
set Sd has more than a single point when for example the geometry is the union of pieces of different 
dimensions [l(i| . In what follows we will look at an upper bound for the dimension: It may happens 
that Tr(|D|~ Q ) = for large enough a, whereas for small values of a, Tr(|D| _a ) = oo. Eventually, 
there is a value of a (say, a c ) for which Tr(|D|~ Qc ) is finite and non-zero. The upper bound for the 
dimension will be S = a c . 



V. COMPUTING (6) 

In order to estimate (S) by means of l|10|) , we rewrite (0 and (J5J as integrals over the eigenvalues 



Xk of D. The procedure is well-known |||, and leads to (C n = ir ' 2 / J|^ =1 kl) 

1*00 

zn{(3) = C n / K l A]A 2 (A fc )e-^"=i A - = C n * ni/J , (11) 



(O(Xi)) n ,0 = J Jd"A]0(A t )| ^ n /^, A 2 (A fc )e-^^^| (12) 

[d n \]0{X t )VnAXk) , (13) 
where A(Afe) = Il < (^i — ^i) ^ s the Vandermonde determinant (Jastrow factor), and [d n X] = 

In random matrix theory, Vl'n^^) is interpreted as the positional partition function of an en- 
semble of equal charged particles (with positions given by Ai) in two dimensions, moving along an 
infinite line, in thermodynamic equilibrium at temperature 7 - the so-called "Dyson gas" |9j. Then, 
"Pn,f)(Xi, A n ) defined in i|13|) is the probability of finding one particle at Ai, one at A2, etc. The 
value of ^>n,p(l) is known from the Selberg's integral. 

In the region (3 < (3 C the partition function Z(f3) is dominated by 00— dimensional matrices. Thus, 
one may try to select the 00— dimensional matrices out of the whole ensemble, and then compute the 
mean of the £-functions following (|10|) . However, from the standpoint of our statistical approach this 
procedure does not seems natural since the sum over n is a key ingredient in the whole construction. 
Hence, we look for a quantity related to the ^-function that captures the statistical nature of our 
model. Let us compute the mean value 

(Tr K \D\-«) n> = /j2 \X k \- a 9(\X k \ k)\ . (14) 

\fc=l / n ,p 



After some considerations one arrives at 



2 ,„ oS f 1 dy 



(Tr Kn JD\- a ) n ^ « -(2n) x -f/3t / JL ^T^. (15) 



G 



The asymptotic behavior of (Tr Kri _ \D\~ a ) n ^ does not depend in an essential way on the particular 
choice of e, as long as we keep e ^ 0. 

Now we use the asymptotic formula (|15|) and search for the value a c for which, as N — > oo and 
(3 — » f} c , (Tr Kn fi \D\~~ a ) diverges (converges to zero) if a < a c (a > a c ), with (Ti Kn P \D\~ ac ) finite 
and non-zero. This gives an upper bound for the dimension of the "condensed" manifold in the 
infinite phase ( (3 < (3 C ), which is (6) < a c . We obtain: 

JV 

(Tr Kn>0 \D\- a )(f3) ~ lim ]T P(n,/3)n x -^ . (16) 

n=l 

In the finite phase ((3 > /3 C ) the sum in Ijltifl converges for a > 0. We conclude that a c = (i.e. 
(S) = 0) for j3 > (3 C , as expected. From the behavior of P(n, /3) in the infinite phase it follows 
that the convergence of the sum in (|16|) is dictated by the behavior of T n a = n 1 ~ : i in the limit 
n N — > oo. For (3 < j3 c we get T n:a — > oo if a < 2, and r n>a — » if a > 2. For a = 2 it turns out 
that (Tr Kn . |Z)| _2 )(/3) ~ 1. Therefore, we obtain the upper bound (5) < 2. 



VI. ONE DIMENSIONAL MODEL 



We would like to apply the ideas of last sections to spaces of dimension one instead of dimension 
zero. Let us consider a collection X of n one dimensional intervals. The corresponding spectral 
triple will be {Ax, Hx, D) where Ax is the algebra of continuous functions on X and Hx — L 2 {X). 
The analogue of the Dirac operator will be the momentum operator —i-§^- We will keep X fixed 
and fluctuate D. 

Let us consider a simple example where X is a pair of disjoint intervals The intervals will 

be parametrized by coordinate x £ [0, 2n]. An element ip £ Hx is a pair of functions ipi(x), ip2(x), 
ipi : Ii — > C and the scalar product is 

C0,x) = / dx(ip$ X i (17) 
Jo 

We can write ip as a column vector and the operator D as the following matrix 

We have not fixed completely the spectral triple. The operator D is fixed only up to boundary 
conditions (BC's) or self-adjoint extensions. The most general BC can be written as 



V>i(2tt)\ = ( 9n 9i2 \ (MO) 

02(27!") J \ 521 522 / V ^2(0) 



(19) 



where 9 is a matrix in U(2). 

We may ask what geometrical properties of X is determined by such BC's. The point of view 
taken by Balachandran at all in |5| is that a BC fixes the global topology. A couple of examples 
will illustrate their point of view. First let us suppose gij = Sij . In this case ipi (0) = ipi (27r) . The 
two intervals fold into a pair of independent circles. Now suppose 511 = 522 = and 1712 = .921 = 1- 
Therefore -0i(27r) = ^2(0) and 02(271") = 0i(O). As one can see, the two intervals are connected to 
make a single circle of size Air. For a generic BC, however, topology is not localized as in these two 
examples but it is rather a superposition of both. We refer to for further details. 
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We would like to compute the partition function 

Z x (P)=J2^ TrD2 (20) 

D 

and look at the probability distribution for the BC's. For finite values of /3 the topology, as given 
by the BC, will fluctuate. We have evidence, however, that the BC gets localized as we increase /3 
and becomes a pair of circles in the infinite limit. The significance of this fact and generalizations 
of this one dimensional model will appear in Q . 



VII. FINAL REMARKS 



We proposed a discrete model for Euclidean quantum gravity, i.e, random geometry, based on 
the framework of noncommutative geometry. The first model contains the mean number of points, 
(n), and the dimension of the space-time, (S), as dynamical observables. We have shown that the 
discrete model has two phases: a finite phase with a finite value of (n) and (5) = 0, and an infinite 
phase with a diverging {n) and a finite (S) ^ 0. An upper bound for the order parameter (5) was 
found, (S) < 2. We also considered the simplest example of another model where dimension is fixed 
but topology fluctuates. In the limit of infinite /3, however, topology gets localized . 
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